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^h Abstract. Photoproduction of p° mesons in ultra-peripheral Pb+Pb collisions has been 

\f*} studied by the ALICE Collaboration at the CERN LHC. The strong photon flux associated 

£S| with relativistic charged nuclei leads to a very large cross section for exclusive photoproduction 

of p° meson in interactions of the type Pb + Pb — > Pb + Pb + p° . For a p° produced at mid- 
rapidity at the LHC, the photon- nucleus center of mass energy is higher than in any previous 
experiment. 

The ALICE detector is a general purpose detector dedicated to study heavy-ion collisions. 
*\ ALICE has excellent performance in the low pr region, and can reconstruct charged particle 

^ tracks with 0.1 GeV/c < pr < 100 GeV/c. In this analysis all tracks were required to be within 

£^ ALICE'S central barrel. Analysis of data from the first heavy ion run at the LHC in 2010 will 

be discussed in this paper. 

> 

o 

00 1. Introduction to ultra— peripheral collisions 

Ultra-peripheral collisions are collisions between hadrons, they can be protons or nuclei, where 
they geometrically miss each other. This implies that the impact parameter is larger than two 
times the radii of the colliding hadrons, see Figure 1. For light vector meson production in 
Pb-Pb collisions at LHC energy, the median impact parameter is usually in the range 200 - 
300 fm. Because such impact parameters are much larger than the short range of the strong 
interaction, the interactions will be mediated by the electromagnetic field. The electromagnetic 
field of a moving charged particle can be treated as a flux of virtual photons. This model is 
referred to as the Weizacker- Williams method. The intensity of the electromagnetic field, and 
therefore the number of virtual photons in the cloud is proportional to Z 2 , where Z is the charge 
of the particle [1]. 

One can divide ultra-peripheral collisions into two categories: two-photon and photonuclear 
interactions. The photonuclear interactions can be further subdivided into coherent and 
incoherent interactions. In coherent photonuclear interactions, the photon interacts with the 
whole target nucleus coherently. In most cases both nuclei will remain intact after the interaction. 
For photonuclear interactions a photon from one nucleus is interacting with the other nucleus 
and, because the photon has spin J p = 1~, a vector meson can be produced. The transverse 
momentum of the resulting particles will be determined by the nuclear form factor, the average 
transverse momentum will be (pt) ~ 60 MeV/c. 

In an incoherent interaction a photon from one nucleus interacts with a single nucleon in the 
target nucleus. This will cause the target nucleus to break up. The transverse momentum is 
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Figure 1. The impact parameter, 6, is larger 
than two times the radius, R, of the colliding 
hadrons. The photon flux from each nucleus 
is proportional to Z 2 . 



higher in incoherent interactions than in coherent interactions, (pt) ~ 400 MeV/c for the p° [6]. 

In two-photon interactions one photon from each of the two colliding nuclei collides and 
produces a fermion pair. An example is given by Pb + Pb —¥ Pb + Pb + e + e _ , where the 
interaction is 7 + 7 — >► e + e~. 

For vector meson production in photonuclear coherent interactions the event will contain 
only two tracks, and the detector will otherwise be empty. There will be one positive and one 
negative track, and the transverse momentum of the track pair will be ^2pt J$ 100 MeV/c [6, 9]. 



2. Model predictions 

The center-of-mass energy per nucleon in Pb-Pb collisions at the LHC in the 2010 and 2011 
Pb-Pb runs was ^snn = 2.76 TeV. This corresponds to a Lorentz factor of 7^ = 1470 for each 
beam in the center of mass system. The p° is a broad resonance with mass M p — 775 MeV/c 2 
and width Y p — 149 MeV/c 2 . At mid-rapidity, y = 0, the mass 775 MeV/c 2 corresponds to a 
mean photon energy of E 1 — 775 MeV/2 = 338 MeV. The equivalent 7-proton center-of-mass 
energy is then 
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y/4:E^m p j L = 43.2 GeV. 
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This is three times higher than at RHIC, and higher than in any fixed-target experiment. 

Three different models predict the cross section for p° production at LHC energies. The model 
by Frankfurt, Strikman, and Zhalov (GGM (Gribov-Glauber Model)) [3, 4] uses a generalized 
vector dominance model in the Gribov-Glauber approach. It includes non-diagonal transitions, 
where the photon fluctuates to a //, but appear as a p° after scattering off the target nucleus. 
The cross section a(p + nucleon) from the Donnachie-Landshoff model, which is in agreement 
with HERA and lower energy data, is used as input for the photonuclear calculation. 

The model developed by Goncalves and Machado (GM) [5] is based on the color dipole 
model in combination with saturation from a Color Glass Condensate approach. Starlight [6, 7] 
is a Monte Carlo event generator developed by Klein and Nystrand. Experimental data for 
7 + p —> p° + p is used in combination with a Glauber model, neglecting the elastic scattering. 
The predictions from the different models can be seen in Figure 2. 
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Figure 2. Model predictions for p° 

photoproduction at central rapidities (y = 0), 
for the models GGM [3, 4], shown with a full 

( ) red line, GM [5], shown with a yellow 

chained ( — • — ) line and Starlight [6, 7], 

shown with a blue dashed ( ) line. (Color 

online.) 



3. The ALICE detector 

The ALICE detector (Figure 3) is a general purpose detector at the CERN LHC [8]. Its main 
goal is to study ultra-relativistic heavy-ion collisions. It consists of a central barrel, a forward 
muon arm, and some other smaller forward detectors. The central barrel has a pseudorapidity 
acceptance of \r]\ < 0.9 and a transverse momentum acceptance of pt > 100 MeV/c. 




Figure 3. The ALICE detector 



For this analysis mainly the central barrel is used. The Inner Tracking System (ITS) and 
the Time Projection Chamber (TPC) are used for tracking and particle identification. The ITS 
consists of six layers of silicon detectors; the two innermost layers are silicon pixel detectors 
(SPD), the two next are silicon drift detectors (SDD), and the two outermost layers are silicon 
strip detectors (SSD). The ALICE TPC, a cyllindrical gaseous detector with a diameter of 
500 cm and a length of 510 cm, is the main tracking device in ALICE. The 557,568 readout pads 
can provide up to 159 ionization samples for track reconstruction, which can be used to calculate 
the energy loss (dE/dx) of the track. The dE/dx of a track is calculated as the truncated mean 
of the dE/dx of the clusters associated with the track. The truncated mean is used in order 
to reduce the fluctuations in cluster energies resulting from the Landau tail [13]. For triggering 
the SPD, the Time-of-Flight detector (TOF) and the VZERO counters are used. The TOF 
detector, which surrounds the TPC, is composed of multigap resistive plate chambers which 
provide an intrinsic resolution of approximately 80 ps. The VZERO counters, on each side of 



the interaction point, consists of 32 tiles of scintilators. They are in this analysis used to define 
rapidity gaps. The VZERO-A has an acceptance of 2.8 < 77 < 5.1, and the VZERO-C has an 
acceptance of —3.7 < 77 < —1.7. To count neutrons from nuclear break up, the Zero Degree 
Calorimeters (ZDC), which are hadronic calorimeters, located ±116 meters on each side of the 
interaction point, are used. 

4. Analysis of 2010 Pb-Pb data 

4.1. Data sets and cuts 

For this analysis data recorded during the 2010 Pb-Pb run, at center-of mass energy ^/snn = 
2.76 TeV, are used. In the first part of the run, when the luminosity was low, a trigger which 
required at least two hits in the TOF detector was used. Later in the run, when the luminosity 
was higher, a tighter trigger definition was needed. A trigger which required at least two hits 
in the TOF detector, at least two hits in the Si-pixel detectors and no activity in the VZERO 
counters was implemented for the last part of the run. 
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Figure 4. Energy loss in the TPC plotted 
versus the momentum of the track, for p° 
candidates 



Figure 5. Energy loss of the positively 
charged particles vs. the energy loss of the 
negatively charged particles in the TPC, for 
p° candidates. 



Based on the characheristics of UPC events, a set of cuts are applied to the data at analysis 
level. The events are required to satisfy one of the two UPC triggers at hardware level. The event 
must have a primary vertex, where the z-position (along the beam axis) is within 10 cm from 
the center. There must be exactly two accepted ITS+TPC tracks, and the VZERO detectors 
should be empty. The two tracks are required to produce ionization in the TPC consistent with 
pions (Figure 4 and 5). The rapidity of the mother particle should be \y pa ir\ < 0.5 and the 
transverse momentum required to be below p^ lT < 150 MeV/c, to get the coherent events. In 
the end it was required that the tracks have opposite charge, and the like sign background (less 
than 2%) was subtracted. The cuts, and the number of events surviving each of them, are listed 
in Table 1 and Table 2. 



4-2. Ace x Eff correction 

To correct for acceptance and efficiency a flat invariant mass 7r + 7r - simulation is used. The 
flat simulation is chosen because the shape of the p° peak would give low statistics for the 
correction at the tails of the distribution. It is assumed that the 7r + 7r - pairs are emitted from 
a transversely polarized parent, as expected for coherently produced p°'s. The generated 7r + 7r - 
events are processed by the ALICE simulation and reconstruction framework using the Geant 



Table 1. Number of events surviving the cuts 
for the TOF-trigger. 



Table 2. Number of events surviving the cuts 
for the SPD+TOF+VZERO vetoed-trigger. 



Cut 


Events left 


Cut 


Events left 


Triggered events 


1,332,041 


Triggered events 


121,487 


Primary vertex 


850,409 


Primary vertex 


103,480 


Two accepted tracks 


47,978 


Two accepted tracks 


26,217 


\V Z \ < 10 cm 


43,413 


\V Z \ < 10 cm 


24,020 


VZERO veto 


8,848 


VZERO veto 


17,567 


PID cut 


7,588 


PID cut 


15,377 


\y\ <o.5 


5,887 


\y\ <0.5 


11,928 


f T atr < 150 MeV/c 


2,749 


f™ r < 150 MeV/c 


6,195 


Unlike sign pairs 


2,699 


Unlike sign pairs 


6,101 


Like sign pairs 


50 


Like sign pairs 


94 



transport to simulate the detector response. These simulated events are then passed through 
the same analysis as is used on the data. The (Ace x Eff) is defined as the ratio of the number 
of reconstructed to generated selected events in the rapidity interval \y\ < 0.5 and pt < 150 
MeV/c as a function of invariant mass in the range 2/77^ < M{ nv < 1.5 GeV/c 2 . 

4-3. Fitting the invariant mass distribution 

The invariant mass distribution is corrected for (Ace x Eff) and fitted with a Breit-Wigner 

function with continuum correction (Equation 2). 
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where 

T(M^) = r • (M p o/M™) x [(M^ - 4ml)/(M 2 p0 

is the momentum dependent width of the p meson, A is the amplitude of the Breit-Wigner 
function and B is the amplitude of the direct non-resonant tt + tt~ production. This function has 
previously been used by the STAR [9] and HI [10] collaborations. 

The fitted invariant mass distribution is shown in linear scale in the left panel of Figure 6, 
and in logaritmic scale in the right panel. The mass of M p o = 767.8 ±3.5 MeV/c 2 and the width 
of T p o = 154.1 ±8.7 MeV/c 2 , are compatible with the PDG values of M = 775.49±0.34 MeV/c 2 
for the mass, and T = 149.1 ± 0.8 MeV/c 2 for the width. 



4-4- Subtraction of incoherent contribution 

The transverse momentum cut pt < 150 MeV/c will leave mostly the coherent events, but 
also some incoherent events will remain. To account for this, one has to find the fraction of 
incoherent events with pt < 150 MeV/c. To do this the Starlight [6, 7] particle generator is 
used. The detector response is also simulated, and the simulated particles are reconstructed 
with the ALICE analysis framework [11]. 

The data points and the simulated distributions are shown in Figure 7. A coherent and an 
incoherent simulated sample was generated. The two samples were scaled to fit the data. The 
contribution from incoherent events under the coherent peak (pt < 150 MeV/c) is found to be 
~ 7%. The coherent peak appears to be slightly narrower in data than in the simulation. This 
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Figure 6. The invariant mass distribution fitted with a Breit-Wigner function with continuum 
correction. Linear (left) and logartimic (right) scale. The data points are marked with a full 

circle (• ), the blue dashed line ( ) is the resonance plus the continuum, the full green line 

( ) is the resonance plus the continuum in the range used to make the fit and the red dotted 

line ( ) is only the resonance. Statistical errors are shown. (Color online.) 
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Figure 7. The transverse momentum distribution compared to Starlight predictions for 
coherent and incoherent photoproduced p°'s. The data points are marked with a full circle 

(• ), the simulated coherent production is marked with a red full line ( ), the simulated 

incoherent production is marked with a full blue line ( ) , and the sum of the incoherent and 



coherent simulation is marked with a dashed green line (- 



(Color online.) 



could be a detector effect or related to the implementation of the nuclear form factor in the 
Monte Carlo. 



4.5. Nuclear break up 

Exchange of additional photons may lead to coherent vector meson production in coincidence 
with nuclear break up [12]. To count the number of neutrons emitted from the nuclei, the 
ZDC neutron counters, located at 116 meters on each side of the interaction point, are used. 



Figure 8 shows the energy deposited in each of the ZDCs, for triggered events with two accepted 
tracks, and the z-position of the primary vertex within 10 cm from the center. The peaks 
corresponding to zero and one neutron in the ZDCs can be clearly identified. The peak centered 
around E — 1380 GeV, which is the beam energy per nucleon, represents one neutron detected. 
In the final selection one distinguishes between the two cases of zero neutrons and one or 
more neutrons in the ZDC. The separation between no neutrons and one or more neutrons is 
defined from the minimum between the first two peaks in Figure 8. As expected, the energy at 
this minimum corresponds to half the beam energy per nucleon. This distinction will make it 
possible to meassure the cross section for photoproduction of p° with and without nuclear break 
up seperately, and this can be compared with model predictions [12]. 
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Figure 8. Energy deposited in the ZDC on the A-side (left) and the C-side (right). Each peak 
is fitted with a Gaussian. The first peak represent zero neutrons, the second peak one neutron 
and so on. Because a pedestal value is subtracted from the signal, the energy goes below zero. 
The data points are marked with a full circle (• ), the fit to each peak each marked with a green 
dashed line ( ), and the sum of the fits is marked with a red full line ( ). (Color online.) 



5. Conclusions and outlook 

The exclusive photonuclear production of p° has been studied by the ALICE collaboration. The 
analysis cuts developed for ultra-peripheral collisions allow a clear separation of the signal, with 
a background contribution estimated from the like-sign contribution of less than 2% (cf. Tables 
1 and 2). The invariant mass and transverse momentum distributions generally agree with the 
expectations, although the coherent peak in the pt distribution appears slightly narrower in 
data. The measurement of the cross section at mid-rapidity will allow a better understanding 
of p° photoproduction. 
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